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Abstract: We study the properties of the space-time that emerges dynamically from the 
matrix model for type IIB superstrings in ten dimensions. We calculate the free energy and 
the extent of space-time using the Gaussian expansion method up to the third order. Unlike 
previous works, we study the SO(d) symmetric vacua with all possible values of d within 
the range 2 < d < 7, and observe clear indication of plateaus in the parameter space of the 
Gaussian action, which is crucial for the results to be reliable. The obtained results indeed 
exhibit systematic dependence on d, which turns out to be surprisingly similar to what 
was observed recently in an analogous work on the six-dimensional version of the model. 
In particular, we find the following properties: i) the extent in the shrunken directions is 
given by a constant, which does not depend on d; ii) the ten-dimensional volume of the 
Euclidean space-time is given by a constant, which does not depend on d except for d = 2; 
iii) The free energy takes the minimum value at d = 3. Intuitive understanding of these 
results is given by using the low-energy effective theory and some Monte Carlo results. 

Keywords: Matrix Models, Superstring Vacua. 



Contents 

1. Introduction 1 

2. The model and the order parameter 4 

3. The Gaussian expansion method 5 

4. Results 9 

5. Summary and discussions 16 

A. Details of the ansatz 17 

B. Free energy from the Krauth-Nicolai-Staudacher conjecture 31 



1. Introduction 

In recent years a lot of efforts have been devoted to constructing a stable (or long-lived 
meta-stable) vacuum in string theory, which is appropriate to describe our real world. 
In this kind of approach, however, there is no objective criterion to pick up one of the 
tremendously many possible vacua, which is commonly referred to as the landscape problem 
in the literature. It would be of course scientifically more desirable if a unique vacuum is 
chosen nonperturbatively and the chosen vacuum actually describes our real world. The 
aim of the present paper is to show that this possibility should not be totally forgotten as 
theorists' dream. Our explicit calculations rather suggest that we might be quite close to 
it although a big twist may still be needed to achieve the final goal. 

In order to address such an issue, we certainly need a nonperturbative formulation of 
superstring theory. There are actually quite a few proposals made in the late 90s [1, 2, 3] 
after the discovery of D-branes. In this paper we study the type IIB matrix model [2], 
which is proposed as a nonperturbative formulation of type IIB superstring theory in 
ten dimensions. This model can be obtained formally by taking the zero- volume limit 
of SU(iV) super Yang-Mills theory in ten dimensions. 1 The integer N can be viewed as 
a sort of regularization parameter, which should be taken to infinity eventually. In this 
formulation of superstring theory, the ten-dimensional target space is represented by the 
ten bosonic Hermitian matrices [4], which originate from the ten-dimensional gauge field 

1 This relationship to SU(iV) super Yang-Mills theory is considered at the level of classical action. There- 
fore the well-known gauge anomaly in the ten-dimensional super Yang-Mills theory is not an issue. 
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in the super Yang-Mills theory. The model has manifest SO(10) invariance, 2 which is 
inherited from the super Yang-Mills theory. The question we would like to ask is whether 
this SO(10) symmetry is spontaneously broken in the large- N limit down to SO(4), which 
is the symmetry corresponding to our space-time. 3 

In ref. [6] two of the authors applied the Gaussian expansion method to this issue. The 
free energy of the SO(d) symmetric vacua for d = 2,4,6,7 was calculated up to the third 
order of the expansion, and d = 4 was found to give the smallest value. Moreover, the 
extent of space-time in the d directions was found to be larger than those in the remaining 
(10 — d) directions. This result motivated higher order calculations up to the eighth order 
[7, 8, 9, 10, 11] for the d = 4 and d = 7 cases. While the results revealed interesting 
qualitative differences between the two cases, the situation turned out to be obscure. 

In the Gaussian expansion method, the result depends on free parameters introduced 
in the Gaussian action. The crucial point is that one can still make a reliable prediction by 
identifying the "plateau region" in the parameter space, in which the result becomes almost 
constant. A practical approach is to obtain the points in the parameter space at which 
the result becomes stationary. As one goes to higher orders, one typically obtains more 
and more stationary points giving totally different results. However, if it turns out that 
there are quite a few points that give approximately the same results, one may regard it 
as indication of the plateau region. Such behaviors were indeed observed in various simple 
models and the obtained results confirmed the validity of the method. (For example, see 
refs. [7, 12, 13].) However, in the case of type IIB matrix model, the plateau region has 
not yet been identified unambiguously. 

Recently it was suggested [14] that the above situation might be due to the extra 
symmetry E^, which was imposed on the shrunken (10 — d) directions in order to make the 
calculations feasible. By imposing a symmetry SO(d) xS^C SO (10), which is stronger 
than just SO(<f), one can reduce the number of free parameters considerably. While this 
makes it much easier to obtain stationary points, it also makes the available stationary 
points strongly restricted by the chosen extra symmetry As a result, one might miss 
the opportunity to observe clear indication of plateaus for the SO(d) symmetric vacua that 
might otherwise be there. 

This possibility was noticed in a similar study for the six-dimensional version of the 
type IIB matrix model [14]. The model can be obtained formally by the zero- volume limit 
of SU(iV) super Yang-Mills theory in six dimensions, and it is expected to have properties 
analogous to the type IIB matrix model such as the spontaneous breakdown of the SO (6) 
symmetry. On the technical side, the analysis becomes much easier than in the type IIB 
matrix model, and it was possible to perform calculations for all the values of d within the 
range 2 < d < 5. Moreover, calculations were done imposing only the SO(d) symmetry 
up to order 3 for d = 3,4,5, and up to order 5 for d = 4,5. Quite a few stationary 

2 This should be considered as a virtue of this model, which is important for the issue we address in this 
paper. In any other proposal for a nonperturbative formulation of superstring/M theory, the space-time 
symmetry is not manifestly preserved. 

3 See ref. [5] for discussions on the idea of the emergent gravity, which is deeply related to the above 
scenario. 
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points were found for each d, and clear indication of plateaus was observed. This enabled 
reliable predictions for both the free energy and the extent of space-time, which exhibited 
interesting qualitative features summarized as follows. 

i) The extent of space-time in the shrunken directions is given by a constant (r), which 
is independent of d. (universal "compactification" scale) 

ii) The six-dimensional volume of the Euclidean space-time is given by a constant (v = 
£ 6 ), which is independent of d except for d = 2. (constant volume property) 

iii) The free energy takes the minimum value at d = 3. 

Intuitive understanding for these properties is also given in ref. [14]. The properties i) 
and ii) can be understood by considering the low-energy effective theory, which is given in 
terms of a system similar to the branched polymer [4, 15]. The two dynamical scales r and 
£, which characterize the properties i) and ii), respectively, are reproduced numerically by 
Monte Carlo simulation [16]. The property iii) as well as the anomaly for d = 2 in the 
property ii) can be understood from the properties of the fermion determinant [17]. 

With all these new insights, we redo the calculations for the type IIB matrix model 
and study the SO (a!) symmetric vacua for all values of d within the range 2 < d < 7. 
Unlike in the six-dimensional case, it is difficult to perform calculations imposing only the 
SO(d) symmetry. However, the calculations in the six-dimensional case [14] suggest that 
the stationary points obtained in the plateau region have all the variety of symmetries in 
the shrunken directions. Therefore, we perform calculations in the ten-dimensional case 
imposing SO(d) x C SO(10) with various possible for each d. In fact we exhaust all 
the possible extra symmetries that leave not more than five free parameters in the Gaussian 
action. By combining all the solutions obtained in this way, we were able to observe clear 
indication of plateaus already at the 3rd order calculations as we did in the six-dimensional 
case. This should be contrasted to the situation with the previous calculations for the type 
IIB matrix model, where one particular symmetry was chosen for each d, and hence the 
number of solutions was not enough to clearly identify the plateau. 

The free energy and the extent of space-time obtained for each d in the plateau region 
indeed exhibit systematic d dependence analogous to i)-iii) observed for the six-dimensional 
case. This is understandable theoretically given that the low-energy effective theory is de- 
scribed by a similar branched-polymer-like system [4, 15] and that the fermion determinant 
has similar properties [17]. Let us emphasize, however, that the Gaussian expansion method 
knows neither of these facts, and yet it revealed the same qualitative behaviors of the two 
models. Moreover, the free energy obtained for the SO(d) symmetric vacua is actually 
quite close to the value obtained from the formula for the partition function conjectured 
by Krauth, Nicolai and Staudacher [18], for both the 6d and Wd models. We consider these 
as strong evidence for the validity of the present calculations. 

The rest of this article is organized as follows. In section 2 we define the model and the 
observable which serves as an order parameter for the spontaneous breaking of rotational 
SO (10) symmetry. In section 3 we describe the method we use to study the model. In 



-3- 



section 4 we present our results for the free energy and the extent of space-time in the 
SO ((f) symmetric vacua (2 < d < 7). Section 5 is devoted to a summary and discussions. 
In appendix A we present the details of the ansatz we use to study each of the SO(d) 
symmetric vacua. In appendix B we derive the value of free energy from the formula for 
the partition function conjectured by Krauth, Nicolai and Staudacher. 

2. The model and the order parameter 

The type IIB matrix model can be obtained formally by the zero-volume limit of D = 10 
SU(iV) pure super Yang-Mills theory, and its partition function is given by 

d^d^e- Sb ~ Sf , (2.1) 
S h = -^Tr[A^,A u } 2 , (2.2) 



S t = -^TrOMCr^LV^]) . (2.3) 



Here (fi = 1, ■ • ■ ,10) are traceless N x N Hermitian matrices, whereas \l/ a (a = 
1, • • • ,16) are traceless N x N matrices with Grassmannian entries. The parameter g 
can be scaled out by appropriate redefinition of the matrices, and hence it is just a scale 
parameter rather than a coupling constant. We therefore set g 2 N = 1 from now on without 
loss of generality. The integration measure for A^ and Vl/ Q is given by 

N 2 -l 10 N 2 -l 16 

a=l [i=Y v a=l a=l 

where A a ^ and ^ are the coefficients in the expansion 

\ = E^Lr 1 A l Ta etc. with respect 
to the SU(iV) generators T a normalized as Tr(T a T b ) = \5 ab . 

The model has an SO(10) symmetry, under which A^ and ^ a transform as a vector 
and a Majorana-Weyl spinor, respectively. The 16 x 16 matrices are the gamma matrices 
after the Weyl projection, and C is the charge conjugation matrix, which satisfies (T^) T = 
CT M Ct and C T = C. 

In order to discuss the spontaneous symmetry breaking (SSB) of SO(10) in the large- N 
limit, we consider the "moment of inertia" tensor [4, 19] 

T^ = lTr(VV) , (2-5) 

which is a 10 x 10 real symmetric tensor. We denote its eigenvalues as Xj (j = 1, • • • , 10) 
with the specific order 

Ai > A 2 > • • • > Aio • (2.6) 

If the SO(10) is not spontaneously broken, the expectation values (Xj) (j = 1, • • • , 10) 
should be all equal in the large- iV limit. Therefore, if we find that they are not equal, it 
implies that the SO(10) symmetry is spontaneously broken. Thus the expectation values 
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(Xj) serve as an order parameter of the SSB. In ref. [17] it was found that the phase of 
the fermion determinant (or Pfaffian, strictly speaking) favors d(> 3)-dimensional config- 
urations, which have Xj (j = d + 1, ■ ■ ■ , 10) much smaller than the others. This suggests 
the possibility that the SO(10) symmetry is broken down to SO(d) with d > 3. Since the 
eigenvalue distribution of represents the extent of space-time in the type IIB matrix 
model [4], the above situation realizes the dynamical compactification to <i-dimensional 
space-time. 

In general one can obtain supersymmetric matrix models by taking the zero-volume 
limit of pure super Yang-Mills theories in D = 3, 4, 6 and 10 dimensions, where the D = 10 
case corresponds to the type IIB matrix model. The convergence of the partition function 
for general D was investigated both numerically [18] and analytically [20]. The D = 3 
model is ill-defined since the partition function is divergent. The D = 4 model has a real 
positive fermion determinant, and Monte Carlo simulation suggested the absence of the 
SSB of rotational symmetry [21]. (See also refs. [22, 23].) The D = 6 model and the 
D = 10 model both have a complex fermion determinant, whose phase is expected to play 
a crucial role [17, 24, 25, 26, 27] in the SSB of SO(D). 

3. The Gaussian expansion method 

Since there are no quadratic terms in the actions (2.2) and (2.3), we cannot perform a 
perturbative expansion in the ordinary sense. Finding the vacuum of this model is therefore 
a problem of solving a strongly coupled system. It is known that a certain class of matrix 
models can be solved exactly by using various large- iV techniques, but the present model 
does not belong to such a category. 4 The use of the Gaussian expansion method in studying 
large- N matrix quantum mechanics has been advocated by Kabat and Lifschytz [28], and 
various black hole physics of the dual geometry has been discussed [29]. Applications to 
simplified versions of the type IIB matrix model were pioneered by refs. [30]. 

The starting point of the Gaussian expansion method is to introduce a Gaussian term 
So and to rewrite the action S = S*b + Sf as 

S = (S + S)-S . (3.1) 

Then we can perform a perturbative expansion regarding the first term (So + S) as the 
"classical action" and the second term (—So) as the "one-loop counter term". The results 
at finite order depend, of course, on the choice of the Gaussian term So, which contains 
many free parameters in general. However, it is known in various examples that there exists 
a region of parameters, in which the results obtained at finite order are almost constant. 
Therefore, if we can identify this "plateau region", we can make concrete predictions. 
It should be emphasized that the method enables us to obtain genuinely nonperturbative 
results, although most of the tasks involved are nothing more than perturbative calculations 
as emphasized in ref. [31]. 

4 Note, however, that the large-iV limit simplifies the calculation in the Gaussian expansion method 
considerably since it allows us to consider only the planar diagrams. 
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There are some cases in which one finds more than one plateau regions in the parameter 
space. In that case, each of them is considered to correspond to a local minimum of the 
effective action, and the plateau which gives the smallest free energy corresponds to the 
true vacuum. These statements have been confirmed explicitly in exactly solvable matrix 
models [13]. 

As the Gaussian action for the present model, we consider the most general one that 
preserves the SU(iV) symmetry. Note, in particular, that we have to allow the Gaussian 
action to break the SO(10) symmetry so that we can study the SSB of SO(10). In practice 
we are going to restrict the parameter space by imposing the SO(d) symmetry with 2 < d < 
7. The plateau region identified for each d corresponds to a local minimum which breaks 
the SO (10) symmetry spontaneously. By comparing the free energy, we can determine 
which local minimum is actually the true vacuum. 

Making use of the SO (10) symmetry of the model, we can always bring the Gaussian 
action into the form 



So = Sob + Sot , (3.2) 

AT 10 

«Sob = yE M M Tr (^) 2 > ( 3 - 3 ) 
N 16 

S ° f = y E AapTriVcPp) , (3-4) 

a,/3=l 

where and A a p are arbitrary parameters. The 16 x 16 complex matrix A a p can be 
expanded in term of the gamma matrices as 

10 

H,v,p=l 

using a 3-form m^ up . We can then rewrite the partition function (2.1) as 

Z = Z (e- (5 - So) ) , (3.6) 
Zo = j dAd^e~ s ° , (3.7) 

where ( • )o is a vacuum expectation value with respect to the partition function Zq. From 
this we find that the free energy F = — log Z can be expand as 

oo 
k=0 

fo = ~ log Z , 



h = ~ E JkTm k+iC k -i{(S h - S ) k - l (S{) 2l ) c q for k>l , (3.8) 



1=0 



where the subscript 'C' in ( • }c,o implies that the connected part is taken. The expansion 
is organized in such a way that it corresponds to the loop expansion regarding the insertion 
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of the 2-point vertex (—So) as a contribution from the one-loop counterterm. Similarly the 
expectation value of an observable O can be evaluated as 



(0) = (0)o + Y J O k , 



k=l 

° k = E ferW k+iC k -i(0 (S h - S ) k - 1 (S { ) 2l )c,o ■ (3.9) 
1=0 ^ + l >- 

In practice we truncate the series expansion at some finite order. Then the free energy 
(3.8) and the observable (3.9) depend on the free parameters and A a p in the Gaussian 
action. We search for the values of parameters, at which the free energy becomes stationary 
by solving the "self-consistency equations" 

4 F = ' *^ F = ' (310) 

and estimate F and (O) at the solutions. As we increase the order of the expansion, the 
number of solutions increases. If we find that there are many solutions close to each other 
in the parameter space which give similar results for the free energy and the observables, 
we may identify the region as a plateau. 

In actual calculation it is convenient to derive the series expansion (3.8) in the following 
way. First we consider the action 

S = S + eS b + y/eS{ , (3.11) 

Z = j dAd^e -5 , (3.12) 

where e is a fictitious expansion parameter. Next we calculate the free energy in the 
e-expansion as 

oo 

F = -\og Z = Y J ^ k h ■ (3.13) 

fc=0 

Each term f k depends on the parameters M M and m^ v \ in the Gaussian action So- Then 
we substitute these parameters as 

M M -»• (1 - e) M M , m^p -> (1 - e) m^ p , (3.14) 

reorganize the series with respect to e, and set e to 1. In this way we reproduce the 
expression (3.8). The action (3.11) is introduced to obtain the ordinary perturbation 
theory for the first term in (3.1), and the final step (3.14) corresponds to taking account of 
the second term in (3.1) as the one-loop counter term. The main task is to obtain the series 
(3.13), which is nothing more than what is required for ordinary perturbation theory. 5 We 
use a similar procedure to obtain the expansion (3.9) for the observables. 



and the partition function 



5 Further simplification is possible by exploiting the fact that the free energy F is related to the two- 
particle irreducible (2PI) free energy through the Legendre transformation [7]. The number of Feynman 
diagrams decreases considerably by the restriction to 2PI diagrams. While we do not use this technique for 
the 3rd order calculations in the present work, it would be crucial in performing higher order calculations. 
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ansatz 


discrete symmetry for shrunken directions 


symbol 


S0(7) 




o 


S0(6) 


ILiA (*7 o o inliA 
^41^ < ,o,y , iu| i ) 


o 






• 




"-*Z^ i ,o | 1 ) 


A 


S0(5) 


.i(b,7,8,9,lU) 


o 






• 




// j a fa 7 q oi i n\ 


A 




Hj^Ip e 7\~\\ X ILi'iiQ. o i n\ 


A 




^joIp. *7Iq o i n\ X ilj'i.i q o i n\ 
^zi^o, / o,y,iuj N ™ j o^o ) yj iu j 


V 




"-'Z^D, f 1 1 J N £\o.\$ \ I) 


▼ 




z ^ d , / i u ) z ^ o , y 1 1 u ) z o , i j* , \ o , y j* | ± , i u j 


o 


S0(4) 


^2(5,611) x ^2(7,811) x ^2(9,1011) 


o 




llj'ii^ p 7\ X oLiKq o i r\\ X aLtiffx. p t\ /q a in\li^ 
o^OjD, i ) o^o.y, iv) z( \ o ,o ., < f ,\o .,y , ±v f \ ± ) 


• 




^4(5,6, 7, 8|1) X ^2(9,10|1) 


A 




iLi^i p. i c o i n\ X lljoI li f; *7 c o i ri'i 
o^o, / -o,y , luj "-'z^ — 1 1 ,D, t ,o,y , iuj 


▲ 




ilj'il c; <; 7\ X /o o in\ X /zjo /" 1 1 q o i n\ 
0(^0,0, i ) o^o.y , iu_) "- j z^ — i ,o,y, luj 


V 


S0(3) 


•^2(4,5 Ij x ^z(o,7|lJ x ^2(8,911; ^ ^a({4,o},{b,7},{8,9}) 


o 




^z(4,5 1UJ x ^z(b,7|10) x ^z(o,y 1(J) 


• 




X/^Q/f/l ^1 fp Ti f o nlN X lLin{ S A r;"l TIN 






/iq / a k 1 1 nN X ajo / fi t 1 1 n\ X /to /■ c n 1 1 nN 
"-^(^4,0 1UJ N z ^z(D, I 1 1U ) N ^Zl^Ojy IU J 


A 




o( \4,5},{b,7},{o,y}) ^ ^ J z( — 1 1 ,10 J 






^3(4,5,6) x ^3(7,8,9) x ^2({4,5,6},{7,8,9} 1 1) x ^2(-|l,10) 


▲ 




R i^A X iLiliT c 1 1 N X ILi^ifcx inliA X iLili Si Q"L /O inl\ 
0^4,0,0 _) z^/,o|!J z^y.lUij z^-[;,Oj-,-[y,lUj'j 


V 




c; i^A X iLiAii q o inliN X ilj^i \i q o i n'l 
o^4,o,D "-"4^ / -o,y , iu| i j "-"z^ — | i ,o,y , iuj 


T 




^2(4,511) x Z 5 (g 7 8 9 io) X Z 2 (_|i ,6,7,8,9,10) 


O 




IljPI^ P 7 Q O inll\ X ILill 1 R fi 7 Q O 1 fl\ 

"-'D^OjD, r ,o,y,iu ij N "-*z^ — o,u,/,o,y,iuj 


+ 




^2(4,5 1) x ^5(6, 7,8,9,10)4,5) 




S0(2) 


All ^ X f\\-\\ X ^9 ('7 S| 1 N X 2^9 fa 1 n| 1 N 

x ^4({3 ; 4},{5,6),{7 ; 8),{9,10}) 


o 




^4(3,4,5,6|1) x ^4(7,8,9,10|1) x ^2({3,4,5,6},{7,8,9,10}) 


• 




^3(3,4,5) x ^3(6, 7, 8) x ^2({3,4,5},{6,7,8} 1 1) 


A 




x ^2(9,10|l) x Z 2 (_|3,4,5,6,7,8) 






Z 3(3,4,5) x Z 3 ( 6; 7 !8 ) X Z 2 ({3,4,5},{ 6 ,7, 8 }|1) 


▲ 




X^2(9,10|l) x Z 2 (_|9,io) 





Table 1: The list of the S0(d) ansatz (d — 2, ••• ,7), which leave us with not more than 5 
parameters in the Gaussian action. The symbols shown in the right- most column are used in figs. 1, 
2, 4 and 5 to represent the solutions to the self-consistency equations. For the S0(3) x ^2{4,5|i) x 
^5(6,7,8,9,io|4,5) ansatz, we find no solutions in the range — 2 < / < 7 displayed in fig. 1, and hence 
we do not assign any symbol. 

Since we are interested in the large- iV limit, we list up Feynman diagrams in the 
double-line notation, and keep only the planar diagrams when evaluating the free energy 
(3.8) and the observable (3.9). 

There are many free parameters in the Gaussian action; i.e., we get 10 from 
and 120 from m^ up . There are as many self-consistency equations as these parameters. 
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Unfortunately it seems impossible to solve them in full generality. However, it is reasonable 
to expect that some subgroup of SO(10) symmetry such as SO(d) with 2 < d < 7 remains 
unbroken. This allows us to impose the corresponding symmetry on the Gaussian action, 
and hence to reduce the number of parameters. 

In order to study the SO(d) symmetric vacuum, we impose the SO(d) symmetry on 
the Gaussian action. For d > 4, this leads to Mi = ■ ■ ■ = and m^ up = unless fi, v and 
p are three different numbers in {d + 1, ■ ■ ■ , 10}. Therefore we obtain (11 — d) parameters 
from M p and io-dC,3 parameters from m^ up . For d = 3, 771423 can also be non-zero, and 
for d = 2, 777i2 M with /i = 3, 4, • • • ,10 can also be non-zero. Thus the total number of 
parameters is 5, 9, 16, 27, 44, 73 for d = 7,6,5,4,3,2, respectively. For d > 8, the SO(d) 
symmetry enforces m^ up = 0, which makes the Gaussian expansion ill-defined. However, 
the SO(a!) symmetric vacua with d > 8 are expected to exist, and they should be realized 
as a plateau with the SO(d) ansatz (d < 7). See ref. [11] for an explicit example of this 
claim for the SO(10) symmetric vacuum. 

It turned out that solving the self-consistency equations (3.10) is possible for the num- 
ber of parameters not more than 5. Except for d = 7, we therefore impose extra discrete 
symmetries (S^) in the shrunken directions a^+i, • • • , xio, where we require that SO(d) x 
is a subgroup of SO(10). Let us introduce the following notations. First, Z p (j 1) ... )i ) repre- 
sents the group of cyclic permutations acting on x^,- ■ ■ ,Xi p . Similarly, ^{{n,^},— ,{i p ,j p }) 
represents the group of cyclic permutations acting on the sets {x^, x^}, • • • , {xi p ,Xj p }. We 
also combine them with a reflection. For instance, the symbol Zp^... t i p \j lt ... j q ) implies that 
we make a reflection (xj 1 , • • • , Xj q ) i-> (— Xj 1 , • • • , — Xj ) together with a cyclic permutation 
of Xi x ■ ■ ■ ,Xi p . The symbol Z 2 (_|j lj ... j ) implies that we make a reflection without any 
cyclic permutations. Note that the permutation is odd if the integer p is even and the 
number of elements in the set to be permuted is odd. In that case, we have to combine it 
with a reflection in an odd number of directions in order to make the imposed symmetry 
an element of SO (10). Table 1 shows the list of the ansatz that leave us with not more 
than 5 parameters. In the previous works [6, 7, 8, 9, 10, 11] one particular extra symmetry 
was chosen for each d so that the number of parameters is reduced to 3. All the extra 
symmetries used in the literature are included in the list of ansatz we study. See appendix 
A for more detail. 

4. Results 

For each ansatz, we first obtain the free energy up to the third order as a function of the 
free parameters in the Gaussian action. By differentiating the free energy with respect to 
the free parameters, we obtain the self-consistency equations, which we solve numerically 
by Mathematica. The free energy evaluated at each solution is plotted in fig. 1 for the 
SO(d) ansatz (2 < d < 7) described in the previous section. More precisely, we actually 
plot "the free energy density" defined as 




where F = — log Z . 



(4.1) 
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Figure 1: The free energy density (4.1) evaluated at the solutions to the self-consistency equations 
at the orders 1 (Top) and 3 (Bottom). Each symbol represents the largest symmetry in Table 1 
that the solution has. The horizontal line represents the value log 8 — | = 1.32944... obtained from 
the KNS conjecture. The data points surrounded by the dashed lines correspond to the "physical 
solutions" , which indicate the plateau region. 
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Figure 2: The zoom-up of fig. 1 (Bottom) for d = 6,5,4, 3 near the "physical solutions", which 
are surrounded by the dashed line. 
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Figure 3: The free energy density averaged over the "physical solutions" for each d is plotted 
against d. We put error bars representing the mean square error when there are more than one 
physical solutions. The horizontal line represents the KNS value / = log 8 — | = 1.32944 . . ., and 
the dotted line connecting the data points is drawn to guide the eye. 
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The horizontal line log 8 — | = 1.32944... represents the result obtained in appendix 
B from the analytic formula of the partition function conjectured by Krauth, Nicolai and 
Staudacher (KNS) [18]. At order 3, the "physical solutions" identified for each d as de- 
scribed below are surrounded by a dashed line. In fig. 2 we zoom up the region containing 
these solutions for d = 6, 5, 4, 3, which gives the free energy density / ~ 1.5, 1.3, 0.7, —1.4, 
respectively. 

In fig. 3 we show the free energy density obtained by averaging over the physical 
solutions for each d at order 3. We put error bars representing the mean square error 
when there are more than one physical solutions. The result decreases monotonically as d 
decreases from 7 to 3, and it becomes much larger for d = 2. Thus, the SO(3) symmetric 
vacuum gives the smallest free energy density. The d-dependence of the free energy density 
is quite analogous to the one observed in the six-dimensional case [14]. There the value of 
the free energy tends to decrease slightly as one goes from order 3 to order 5. Considering 
such artifacts due to truncation, we speculate that the KNS conjecture actually refers to 
the partition function for the SO (10) symmetric vacuum. 

Let us then explain how we identify the "physical solutions" . For that we also refer to 
the results for the extent of space-time. In appendix A we give the explicit form of {T^) 
for each ansatz. Let us note that it is not diagonal in general, and therefore we need to 
diagonalize it in order to obtain the eigenvalues (Xj). Then we find that some solutions 
give smaller values in the directions involved in the preserved SO(d) symmetry than in 
the remaining directions. Such solutions are not shown in the figures, and will not be 
considered in what follows. For the SO(d) ansatz, the d large eigenvalues (1 < j < d) 
are equal due to the imposed SO(d) symmetry, and we denote the corresponding value as 
R 2 . The remaining (10 — d) eigenvalues for each solution turn out to be quite close to each 
other and we denote the mean value as r 2 . In figs. 4 and 5, we plot the values of R 2 and 
r 2 evaluated at the solutions for each ansatz. 

At order 3, we find a set of solutions for 3 < d < 6 giving similar values for the 
free energy density / and the extent of space-time R 2 and r 2 . We consider this as the 
concentration of solutions [7] , which indicates a plateau region in the space of parameters 
explained in section 3. Thus we can pick up the "physical solutions" for 3 < d < 6 without 
much ambiguity. For the SO(2) ansatz, we find two solutions close to each other with 
/ ~ 3.0, R 2 ~ 3.6 and r 2 ~ 0.11, which we identify as the physical solutions. For the 
SO(7) ansatz, we obtain only two solutions at order 3, which are not close to each other. 
From fig. 1 (Bottom), however, it looks reasonable to identify the one with smaller free 
energy density as the physical solution. 

In fig. 6 we plot the result for R 2 and r 2 averaged over all the physical solutions for 
each d. We put error bars representing the mean square error when there are more than one 
physical solutions. We find that r 2 stays almost constant at r 2 = 0.1 ~ 0.15, which seems 
to be universal for all the SO(<i) symmetric vacua with 2 < d < 7. On the other hand, 
the results for R 2 are found to be larger for smaller d. Let us test whether this behavior is 
consistent with the constant- volume property [14], which is given by R d f 10 ~ d « £ 10 , where 
f 2 represents the universal "compactification" scale, which we leave as a free parameter in 
the present analysis. 
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Figure 4: The extent of space-time R 2 in the extended directions evaluated at the solutions to 
the self-consistency equations at the first order (Top) and the third order (Bottom). Each symbol 
represents the largest symmetry in Table 1 that the solution has. The data points surrounded by 
the dashed lines correspond to the "physical solutions" , which indicate the plateau region. 
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Figure 5: The extent of space-time r 2 in the shrunken directions evaluated at the solutions to 
the self-consistency equations at the first order (Top) and the third order (Bottom). Each symbol 
represents the largest symmetry in Table 1 that the solution has. The data points surrounded by 
the dashed lines correspond to the "physical solutions" , which indicate the plateau region. 



- 14 - 



4 




1 * * 1 1 1 1 1 

2 3 4 5 6 7 

d 

Figure 6: The extent of space-time R 2 and r 2 in the extended and shrunken directions, respectively, 
are plotted against d after taking the average over the "physical solutions" for each d. We put error 
bars representing the mean square error when there are more than one physical solutions. The solid 
and dashed lines connecting the data points are drawn to guide the eye. 
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Figure 7: The extent of space-time R 2 in the extended directions shown in fig. 6 is plotted in 
the log scale against l/d. The straight dashed line represents a fit to the behavior (4.2) excluding 
d = 2, from which we obtain f 2 = 0.155 and H 2 = 0.383. 
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From the constant-volume property, we obtain 

10 f£ 2 \ 

logi? 2 « —log/ -J +logf 2 . (4.2) 

In fig. 7 we therefore plot R 2 in the log scale against 1/d. Indeed we find that the results 
can be fitted to a straight line except for d = 2. From the fit, we obtain £ 2 = 0.383 and 
f 2 = 0.155. Note that the obtained value for f 2 is consistent with r 2 obtained directly 
by the Gaussian expansion method as the extent of space-time in the shrunken directions. 
The obtained value for £ 2 is reproduced by Monte Carlo simulation [16]. 



5. Summary and discussions 

In this paper we have applied the Gaussian expansion method to the type IIB matrix 
model, which was conjectured to be a nonperturbative formulation of type IIB superstring 
theory in ten dimensions. In particular, we have investigated the dynamical properties 
of the model associated with the SSB of SO (10), which was speculated to realize the 
dynamical compactification. Unlike previous works, we studied the SO(d) symmetric vacua 
for 2 < d < 7 systematically. Moreover, we were able to observe clear indication of plateaus, 
which is crucial for the method to be reliable. This was made possible by allowing various 
extra symmetries in the shrunken directions as suggested by the success of a similar work 
on the six-dimensional version of the model. 

Indeed our results bear surprising similarity to the results for the 6d case. The free 
energy was found to decrease as d is lowered until one reaches d = 3, and then it becomes 
much larger for d = 2. This implies that the SO(10) symmetry is spontaneously broken 
down to SO (3). The "compactification" scale is universal for all d. The extent of space-time 
obeys the constant-volume property for d > 3, while the volume for d = 2 is a bit smaller. 
All these properties are common to the 6d and lOd cases, and they can be understood 
intuitively as follows. (See ref. [14] for more detailed discussions.) 

First of all, the long-distance effective theory of these models can be obtained by 
integrating out the off-diagonal elements of the matrices perturbatively [4, 15]. One then 
finds that the effective theory for the diagonal elements is described by a branched-polymer- 
like system, where the diagonal elements of the bosonic matrices are identified with the 
coordinates of the vertices in the polymer. This naturally explains the constant-volume 
property since the branched polymer tends to occupy a fixed volume for entropic reasons. 
The anomaly for d = 2 can be understood as a consequence of the fact that the fermion 
determinant disfavors two-dimensional configurations [17]. The volume that appears in the 
description of the constant-volume property can be calculated by Monte Carlo simulation, 
and it agrees quite well with the value suggested by the Gaussian expansion method for 
both 6d and lOd cases [16]. 

The driving force for the SSB comes from the phase of the fermion determinant [17] as 
has been clearly demonstrated in a simplified non-supersymmetric model [25], which was 
studied by both the Gaussian expansion method [26] and a Monte Carlo method [27] , giving 
consistent results. There the free energy was determined by subtle competition between the 
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effect of the phase, which favors smaller d, and the dynamics of the phase-quenched model, 
which favors larger d. In the present super symmetric matrix models, the latter effect is 
suppressed by 1/N [16, 24], and the free energy is determined essentially by the effect of 
the phase. This explains the monotonic decrease of the free energy for smaller d. The 
anomaly for d = 2 can be understood again from the property of the fermion determinant. 

The extent in the shrunken directions is determined by the fluctuation of the off- 
diagonal elements since the diagonal elements cannot fluctuate in the shrunken directions 
due to the strong effect of the phase. At the leading order of the long-distance approxi- 
mation, the dynamics of the off-diagonal elements is given by the Gaussian action, which 
couples weakly to the dynamics of the diagonal elements. This explains the universal 
compactification scale. 

Since our results are obtained at the 3rd order of the Gaussian expansion, it remains 
to be seen whether the main conclusions change qualitatively at higher orders. In fact the 
6d case has been studied both at the 3rd and 5th orders. The results did change slightly, 
but the qualitative features seemed to be robust. It is nonetheless desirable to perform 
the 5th order calculations in the lOd case as well. It would be also important to perform 
Monte Carlo studies analogous to [27] and to see whether the properties suggested by our 
work can be reproduced. 

Our calculations clearly demonstrate the SSB of SO(10), which is itself an interesting 
dynamical property of the type IIB matrix model. However, in view of the scenario for 
dynamical compactification, there are two problems. One is that it seems more likely 
that the minimum of free energy occurs at d = 3 instead of d = 4, at least within the 
approximation of the present work. The other is that the ratio of the extent of space- 
time in the extended directions and that in the shrunken directions is finite for all d, as 
suggested by the constant volume property and the universal compactification scale. We 
consider that the nontrivial dynamics of the type IIB matrix model suggested by the present 
work confirms the validity of the basic idea to use matrices as the microscopic degrees of 
freedom in formulating superstring theory nonperturbatively. The detailed properties of 
the space-time emerging from the model, however, suggest that it still lacks some important 
ingredient in order to be capable of describing our real world. 
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A. Details of the ansatz 

In this section we describe the details of the ansatz, which are listed in Table 1. We exhaust 
all the symmetries SO(d) x with 2 < d < 7, which reduce the number of parameters 
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in the Gaussian action to 5 or less. We also describe the general form of the moment of 
inertia tensor (2.5) for each ansatz. 

SO (7) ansatz 

• SO(7) 

5 parameters (M, Mg, Mg, Mio, rh) 

M M = (M, • • • , M, Mg, Mg, Mio) , "28,9,10 = rh , and zero otherwise. (A.l) 

7 

The moment of inertia tensor takes the form 



/SO (7) part 



\ 



(A.2) 



C x x 
x C x 
\ xx C"7 

Here and henceforth "x" represents a component which is found to vanish up to 
order 3 by explicit calculation, although there is no such symmetry that enforces it. 

By further imposing the symmetry Z 3 ( 8 g jl0 - ) , which enforces Mg = Mg = Mio, on e 
obtains the "SO (7) ansatz" used in refs. [6, 7, 8, 9, 10, 11]. 

SO (6) ansatz 

• SO(6) x Z 4 ( 7j8)9jl o|i) 
3 parameters (M, M, rh) 

M M = (M,--- .M.M.M.M.M; , 

6 

""17,8,9 = ""18,9,10 = "i-7,9,io = "1-7,8,10 = rh , and zero otherwise. (A. 3) 

The moment of inertia tensor takes the form 



/SO(6) part 



(Tuv) = 



\ 



\ 



C a a a 

a C a a 

a a C a 

a a a C / 



(A.4) 



This ansatz is equivalent to the "SO(6) ansatz" used in ref. [6]. 

• SO(6) x Z 3(8i9)10) 

5 parameters (M, M, M, rh, rh') 

M M = (M, • • • , M, M, M', M', M') , 



"^7,8,9 = "J-7,9,10 = -"i-7,8,io = rh , m 8 ,9,io = rh' , and zero otherwise. (A.5) 
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The moment of inertia tensor takes the form 

/SO(6)part 

<7^> = 



V 



C 



\ 



p p p 



C a a 
a C a 
a a C / 



• SO(6) X Za^am 

5 parameters (M, M, M', M", fh) 



M, 



it 



( M, • ■ ■ ,M,M,M,M',M") , 

6 



m 7,9,io = m 8,9,io = "fh , and zero otherwise. 
The moment of inertia tensor takes the form 



/SO(6)part 



(Tuu) = 



C a 
a C 



x x 



X X 



a 



\ 



C" ) 



SO(5) ansatz 



(A.6) 



(A.7) 



(A.8) 



SO(5) x Z 5 ( 6)7j8)9il0) 

4 parameters (M, M, fh, fh') 

M M = ( M, • • • ,M , M,-- - ,M ) , 

5 5 

"1-6,7,8 = "17,8,9 = "1-8,9,10 = "^6,9,10 = »™6,7,10 = "I , 

^6,7,9 = "T-7,8,10 = w-6,8,9 = "1-7,9,10 = »™6,8,io = m' , and zero otherwise. (A.9) 
The moment of inertia tensor takes the form 



<7» 



SO(5)part 


\ 
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a 


a 
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C 
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c) 



(A.10) 



SO(5) x Z 4 ( 6)7j8)9 | 1 ) 
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5 parameters (M, M, M', rh, rh') 

M„ = (M, ■ • ■ , M, M, M, M, M, M') 



"16,7,8 = "17,8,9 = "16,8,9 = "16,7,9 = "1 , 

"i6,7,io = "17,8,10 = "T-8,9,10 = -"16,9,10 = "i' , and zero otherwise. 



The moment of inertia tensor takes the form 

/ SO(5)part 



(A.ll) 




C a a a 

a C a a 

a a C a 

a a a C 



P P P P 



a 



(A.12) 



SO(5) x Z 4 ( 6j7i8j9 | 10 ) 

5 parameters (M, M, M', m, rh') 



M, 



„ = (Mr ■ ,M, M, M, M, M, M') 



"16,7,8 = "17,8,9 = "16,8,9 = "16,7,9 = "1 , 

"i6,7,io = -"17,8,10 = "i8,9,io = "i6,9,io = rh' , and zero otherwise. 
The moment of inertia tensor takes the form 



(A.13) 



\ 



SO(5)part 




\ 
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a 
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a 
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a 
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a 
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C 



















a) 



(A.14) 



SO(5) x Z 2 ( 6j 7|i) x Z 3 ( 8j9)1 o) 

5 parameters (M,M,M',rh,rh') 



Mp = (M, ■ ■ ■ ,M, M, M, AT, M , M ) , 



"16,8,9 = "16,9,10 = -"16,8,10 = "17,8,9 = "17,9,10 = ""17,8,10 = "I , 

m 8,9,io = rh' i an d zero otherwise. (A. 15) 
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The moment of inertia tensor takes the form 

/ SO(5) part 



V 



C a 
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a C 
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P p 
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p P 
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P p 
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a j 



(A.16) 



SO(5) X ^2(6,718,9,10) x ^3(8,9,10) 

5 parameters (M, M , M',m, m') 



= (M, • • , M, M, M, M', M', M') , 

-""16,8,10 = -"17,8,10 = m , 



"16,8,9 — "17,8,9 — "16,9,10 — m 7,9, 10 

~ / 



TO 6,7,8 = m 6,7,9 = m 6,7,io = m ' > an d zero otherwise. 
The moment of inertia tensor takes the form 



(A.17) 
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SO (5) part 
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C 


a 







a 


C 













a p p 










pop 










P P c) 



(A.18) 



• SO(5) X Z 2 ( 6)7 |i) X Z 2 ( 8 ,9|l) 

5 parameters (M,M,M',M",m) 



M M = (M, • • • . M, M. M, M', M', M") , 

5 

"i6,8,io = "i7,8,io = "i6,9,io = ni7,9,io = "i , and zero otherwise. 
The moment of inertia tensor takes the form 



(A.19) 



/ SO (5) part 



C a 
a C 



x x 

X X 



X X 



C p 
p c 



X X 



C" 



(A.20) 



SO(5) X Z 2 ( 6)7 |io) X Z 2 ( 8) 9| 10 ) X Z 2 ({ 6 ,7},{8,9}|1,10) 
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5 parameters (M, M, M', rh, m!) 

M M = (M, • • • ,M, M, M, M, M, M') , 

S v ' 

5 

"26,8,10 = -"17,8,10 = -""16,9,10 = "17,9,10 = fk , 

"i6,7,io = — "18,9,10 = "i' j an d zero otherwise. 
The moment of inertia tensor takes the form 
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a) 



(A.21) 



(A.22) 



SO (4) ansatz 



SO(4) X Z 2 (5, 6 |i) X Z 2 ( 7j8 |i) X Z 2 ( 9 ,io|i) 

5 parameters (M,M,M',M",m) 



M,, = (M, M, M, M, M, M, M', M', M". M") , 

"15,7,9 = "16,7,9 = "15,8,9 = "16,8,9 = "15,7,10 = "16,7,10 

= "15,8,10 = "i6,8,io = "i , and zero otherwise. 



(A.23) 



The moment of inertia tensor takes the form 





/SO (4) part 
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C" 


J 



(A.24) 



By further imposing Z3({5 i6 },{ 7 ,8},{9, 10 }), which enforces M = Ad' = M", one obtains 
the "SO (4) ansatz" used in ref. [6]. 



SO(4) X Z 3 ( 5 , 6)7 ) X Z 3 ( 8 ,9, 10 ) X Z 2 ({ 5 , 6 ,7},{ 8 ,9,io}|l) 
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4 parameters (M, M, rh, rh') 

M„ i .\/. M, M, M, M, ■ ■ ■ ,M) , 



"15,6,7 = "18,9,10 = "1 , 
"^5,6,8 = "15,6,9 = "15,6,10 = "16,7,8 = "^6,7,9 = "16,7,10 
-Wl5,7,8 = ""15,7,9 = -"15,7,10 = "15,8,9 = "16,8,9 = "17,8,9 



— "15,9,10 — "16,9,10 — "17,9,10 — ~ "15,8,10 — —"16,8,10 

and zero otherwise. 
The moment of inertia tensor takes the form 



""17,8,10 



m 



(A.25) 



/ SO (4) part 




\ 




C a a 


P P P 




a C a 


P P P 




a a C 


P P P 




P P P 


C a a 




P P p 


a C a 


V 


P P P 


a a C J 



(A.26) 



By imposing SO(3) x SO(3) instead of Z3 x Z3, which enforces rh' = 0, one obtains 
the "SO(4) ansatz" used in refs. [7, 8, 9, 10, 11]. 

• SO (4) X ^ 4 (5,6,7,8|l) x Z 2(9,10|l) 

5 parameters (M, M, M', rh, rh') 

M„ = (M, M, M, M, M, M, M, M, M' , M') , 

"15,6,9 = "16,7,9 = "17,8,9 : 

"15,6,7 = "16,7,8 = "15,7,8 = "15,6,8 = m 
The moment of inertia tensor takes the form 



"15,8,9 = "15,6,10 = "16,7,10 = "17,8,10 = -"15,8,10 
and zero otherwise. 



= m , 
(A.27) 



/SO(4) part 



<7» 



C a a a 


P 


P 


a C a a 
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a a C a 
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P P P P 
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P P P P 
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C'J 



\ 



V 

SO(4) x Z 5 ( 6)7j8)9il0 ) x Z 2 (_|i j6) 7 j8)9j io) 
5 parameters (M, M, M', rh, rh') 

M„ = (M, M, M, M, M,M', ■ ■ ■ , M') 



(A.28) 



"15,6,7 = "15,7,8 = "15,8,9 = "15,9,10 = -"15,6,10 = m , 
"15,6,8 = "15,7,9 = "15,8,10 = -"15,6,9 = -"15,7,10 = rh' , and zero otherwise. (A.29) 
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The moment of inertia tensor takes the form 



/SO(4) part 



<7» 



C 



\ 



X X X X X 



C a a a a 

a C a a a 

a a C a a 

a a a C a 

a a a a C J 



SO(4) 



3(5,6,7) x ^3(8,9,10) x ^2(-|l, 8,9,10) 



5 parameters (M, M , M',rh,fh') 



M,, = (M, M, M, M, M, M, M, M', M' , M') , 

"15,8,9 = "15,9,10 = — "15,8,10 = m 6,8,9 = "16,9,10 = — "16,8,10 
= "1-7,8,9 = "1-7,9,10 = -"1-7,8,10 = rfl , 
and zero otherwise. 



"15,6,7 = rhf , 



The moment of inertia tensor takes the form 

/SO (4) part 



(Tau) = 
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C 13 13 











(3 C 13 











P P C'J 



SO(3) ansatz 



SO(3) X Z 2 (4 )5 |i) X Z 2 ( 6j7 | 1 ) X Z 2 ( 8)9 |i) X Z 3 ({4 j 5} j { 6j7 } j { 8j9 }) 

5 parameters (M, M ', M' ,rh,rh') 

M„ !.U..U..\/. .!/.••• ,M,M') , 



"14,6,8 = "15,6,8 = "14,7,8 = "^5,7,8 = "14,6,9 = »™5,6,9 = "14,7,9 = "^5,7,9 = 
"14,6,10 = "15,6,10 = "14,7,10 = "15,7,10 = ""14,8,10 = —"15,8,10 
= -"14,9,10 = -"15,9,10 = "16,8,10 = "17,8,10 = "16,9,10 = "17,9,10 = ™! 

and zero otherwise. 



The moment of inertia tensor takes the form 

/SO(3) part 



V 



C a 
a C 
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(3 13 
f3 f3 
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C a 
a C 



7 7 



(A.34) 



C'J 



• SO(3) X Z 2 (4 ) 5|i ) X Z 2 ( 6 ,7|10) X Z 2 (8,9|10) X %3({4,5},{6,7},{8,9}) x Z 2 ({4,5},{6,7}) 

5 parameters (M, M, M', m, rh) 

M, i M. M. M. M . ■ ■ ■ .M..M') , 

6 

"11,2,3 = m , ""14,5,10 = "i6,7,io = "i8,9,io = "* , and zero otherwise. (A.35) 
The moment of inertia tensor takes the form 



/SO(3) part 
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(A.36) 



C'J 



• SO(3) X Z 2 (4 )5 |io) X Z 2 ( 6 ,7|io) x Z 2 ( 8)9 |io) X Z 3 ({ 4)5 } i { 6i7 } i { 8i9 }) X Z 2 (_| 1;10 ) 

4 parameters (M, M,M',m) 

M„ !.U..U..\/. .!/.••• ,M,M') , 



"14,6,8 = "15,6,8 = "14,7,8 = "15,7,8 = "14,6,9 = "15,6,9 = "14,7,9 = "15,7,9 = m , 

and zero otherwise. (A. 37) 

The moment of inertia tensor takes the form 

/SO(3) part 



C a 
a C 



x x 

X X 







X X 
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(A.38) 
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S0(3) X ^3(4,5,6) x Z 3(7,8,9) x Z 2({4,5,6},{7,8,9}|1) x Z 2(-|l,10) 

5 parameters (M, M, M', m, fh!) 

M„ i.U..U..\/. .!/.••• ,M,M') , 

6 

m^ 5jV = m 5)6yU = -m 4) 6,i, = m 7 , 8# = "i8,9,/* = -m 7t g u 
where (i = 4, 5, 6 and ^ = 7, 8, 9 , 
^4,5,6 = "17,8,9 = "i-' > and zero otherwise. 
The moment of inertia tensor takes the form 



/SO (3) part 



V 



C a a 
a C a 
a a C 



P P P 
P P P 
P P P 







p p p 

p 13 P 
P P P 
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C'J 



SO(3) X Z 3 ( 4;5)6 ) X Z 2 ( 7j8 | 1 ) X Z 2 ( 9) i |i) X Z 2 ({7 i8 },{9,10}) 

5 parameters (M, M, M', fh, fh') 

M jt = (M, M, M, M, M, M, M', M', M', M') , 

"2-4,5,7 = "14,5,8 = "1-5,6,7 = "2-5,6,8 = -™4,6,7 = -"2-4,6,8 
= "2-4,5,9 = "1-4,5,10 = "1-5,6,9 = "15,6,10 = —"14,6,9 = ""14,6,10 

"1-4,5,6 = fh' j and zero otherwise. 
The moment of inertia tensor takes the form 



/SO (3) part 



\ 



C a a 
a C a 
a a C 



111 



111 
111 



1 1 
1 1 
1 1 



C (3 

b c 



p P 
p P 
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i i 
i i 
i i 



p P 
p p 



a p 
P cj 



SO(3) x Z 3 ( 4j5i6 ) x Z 4 ( 7i8j9i10 | 1) x Z 2 (_|7 j8)9j io) 
5 parameters (M, M , M' ,fh,fh') 

M,, = (M, M, M, M, M, M, M' , M', M', M') , 

"14,7,8 = "14,8,9 = "14,9,10 = -"14,7,10 = "2-5,7,8 = "15,8,9 
= "15,9,10 = -"15,7,10 = "16,7,8 = "16,8,9 = "16,9,10 = -"16,7,10 

' and zero otherwise. 



The moment of inertia tensor takes the form 



V 

SO(3) X Z 2 ( 4)5 |i) X Z 5 ( 6j7)8i9)10 ) X Z 2 (_|i i6) 7 i8 ,9,10) 

5 parameters (M, M, M ', m, to') 



SO(3) part 
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a 



.1/,, = (M, M, M, M, M, M', • • • ,M') , 



"14,6,7 = "14,7,8 = "14,8,9 = "1-4,9,10 = ""14,6,10 
= ?"5,6,7 = "1-5,7,8 = "15,8,9 = "1-5,9,10 = ""15,6,10 = "l , 

"14,6,8 = "14,7,9 = "14,8,10 = ""14,6,9 = ""14,7,10 
"15,6,8 = "15,7,9 = "15,8,10 = ""15,6,9 = ""15,7,10 = m' 
and zero otherwise. 



The moment of inertia tensor takes the form 

/SO (3) part 



V 

SO(3) X ^6(5,6,7,8,9,10|l) x ^2(-|5,6, 7,8,9,10) 

5 parameters (M, M, M', m, in') 
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a p p p p 








p a p p p 








p p a p p 








p p p a p 








p p p p C'J 



= (M, M, M, M, M', ■ ■ ■ ,M') , 



"14,5,6 = "14,6,7 = "14,7,8 = "14,8,9 = "14,9,10 = ""14,5,10 = m 
"14,5,7 = "14,6,8 = "14,7,9 = "14,8,10 = -"14,5,9 = ""14,6,10 = "l 
and zero otherwise. 



The moment of inertia tensor takes the form 

/SO (3) part 



V 

SO(3) X Z 2 (4 )5 |i) X ^5(6,7,8,9,10|4,5) 

5 parameters (M, M , M',m,rh') 
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M p = (M, M, M, M, M, M', ■ ■ ■ , M') , 

5 

"T-6,7,8 = W-7,8,9 = "1-8,9,10 = "1-6,9,10 = "1-6,7,10 = ™, , 
"16,7,9 = "17,8,10 = "16,8,9 = "17,9,10 = "16,8,10 = "i' , 
and zero otherwise. 



The moment of inertia tensor takes the form 

/SO(3) part 
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(A.48) 



(A.49) 



V 

SO(2) ansatz 

• SO(2) X Z 2 ( 3j 4|i) X Z 2 ( 5)6 |i) X Z 2 ( 7)8 |i) X Z 2 ( 9) i |i) X Z 4 ({3 j 4} j {5 )6 } j { 7j8 } j { 9)10 }) 

3 parameters (M, M, m) 



(A.50) 



M^iM.M^M,--- ,M) , 



"13,5,7 = "13,6,7 = "13,5,8 = "13,6,8 = "1-4,5,7 = ^4,6,7 = "1-4,5,8 = "14,6,8 
= "15,7,9 = "15,8,9 = "15,7,10 = "15,8,10 = "^6,7,9 = "16,8,9 = "16,7,10 = "16,8,10 
= "13,7,9 = "13,8,9 = "13,7,10 = "13,8,10 = "14,7,9 = "14,8,9 = "14,7,10 = "14,8,10 
"13,5,9 = "13,6,9 = "13,5,10 = "13,6,10 = "14,5,9 = "H6,9 = m 4,5,10 = "14,6,10 = fh , 

and zero otherwise. (A. 51) 
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The moment of inertia tensor takes the form 
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f D\J yZi) pdll 
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(A.52) 



This ansatz is equivalent to the "SO(2) ansatz" used in ref. [6]. 

• SO(2) X Z 4 ( 3) 4 j5)6 | 1 ) X ^ 4 (7,8,9,10|l) x Z 2({3,4,5,6},{7,8,9,10}) 

4 parameters (M, M, m, in') 

M^ = (M,M, M,--- ,M) , 

8 

"13,4,5 = "14,5,6 = "13,5,6 = "1-3,4,6 = "^7,8,9 = "18,9,10 = "17,9,10 = "17,8,10 = fh , 

"J-3,4,7 = "14,5,7 = "J-5,6,7 = -"13,6,7 = ^3,4,8 = "14,5,8 = "15,6,8 = ""13,6,8 
= "13,4,9 = "14,5,9 = "15,6,9 = ~ "13,6,9 = "1-3,4,10 = "14,5,10 = "15,6,10 = -"13,6,10 
= "1-3,7,8 = "T-3,8,9 = "13,9,10 = ""13,7,10 = "14,7,8 = "14,8,9 = "1-4,9,10 = -"14,7,10 
= "15,7,8 = "15,8,9 = "15,9,10 = ""15,7,10 = "16,7,8 = "16,8,9 = "16,9,10 = ""16,7,10 = "l' 



and zero otherwise. 
The moment of inertia tensor takes the form 



(A.53) 





/SO (2) part 
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(A.54) 



SO(2) X Z 3 ( 3)4j5 ) X Z 3 ( 6j7;8 ) X Z 2 ({3 J 4 )5 } ) {6,7,8}|1) x Z 2(9,10|l) x Z 3(-|3,4,5,6,7,8) 
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5 parameters (M, M, M', rh, rh') 

S.U..U..U.--- ,M,M',M') , 

S v ' 

6 

"13,4,9 = "14,5,9 = -"13,5,9 = "16,7,9 = ^7,8,9 = -"16,8,9 
= "13,4,10 = "14,5,10 = ""13,5,10 = "1-6,7,10 = "1-7,8,10 = -"1-6,8,10 = T, 
"1-3,6,9 = "14,6,9 = "15,6,9 = "1-3,7,9 = ^4,7,9 = 1"5,7,9 
= "13,8,9 = "14,8,9 = "15,8,9 = "13,6,10 = "14,6,10 = "15,6,10 
= "13,7,10 = "14,7,10 = "15,7,10 = "13,8,10 = "14,8,10 = "15,8,10 = rh' 
and zero otherwise. 

The moment of inertia tensor takes the form 





/SO(2)part 
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P C'J 



SO(2) X Z 3 ( 3j 4 i5 ) X Z 3 ( 6 78 ) X ^ 2 ({3,4,5},-[6,7,8}|l) x Z 2(9,10|l) x Z 2(-|9,10) 

5 parameters (M, M, M', rh, rh') 

= (M, M, M, • • • ,M,M',M') , 



"13,4,6 = "14,5,6 = -"13,5,6 = "13,4,7 = "M,5,7 = _m 3,5,7 
= "13,4,8 = "14,5,8 = -"13,5,8 = »™3,6,7 = m 3J:8 = -m 3A8 
= "14,6,7 = "14,7,8 = -"14,6,8 = "^5,6,7 = "^5,7,8 = ""1-5,6,8 = ™, 

"13,4,5 = "i6,7,8 = m' > and zero otherwise. 
The moment of inertia tensor takes the form 



/SO(2)part 



(Tau) = 
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B. Free energy from the Krauth-Nicolai-Staudacher conjecture 



In this section we derive the value of free energy from the analytic formula for the partition 
function conjectured by Krauth, Nicolai and Staudacher (KNS) [18]. This conjecture was 
obtained by combining earlier analytic works [32, 33] with Monte Carlo results at small N 
[18]. For the D = 10 model, the formula reads 

N(N+1) JV-l 

Skns = |p(Sb + Sf) , (B.2) 

where an d Sf are defined by (2.2) and (2.3) respectively, and the sum runs over all the 
divisors of N. The value of this sum is smaller than X^m=i = IT ^ or arbitrary N, and 
hence it does not contribute to the "free energy density" (4.1) in the large- N limit. 

As one can see from (B.2), the definition of the action Skns differs from ours by 
the factor of 2/N. In order to absorb this factor, we introduce the rescaled variables 
A'^ = (2/N) 1 / 4 and ^f' a = (2//V) 3 / 8 \I/ Q , whose integration measure is given by 

dA'd^'=( — \ dAd^ . (B.3) 

As a result, the partition function (2.1) can be obtained as 

Z = (f f"'"' Z KNS = ( n aV'^j. (B.4) 

^ ' \fc=l / m\N 

From this, we obtain the large- N asymptotics for F = — \ogZ as 

-31ogiV+ flogS- f\ +0 mp\. (B.5) 



N 2 -l b V 4 / V^ 2 , 

The Gaussian expansion method reproduces the first term of (B.5) correctly for any ansatz. 
Substituting (B.5) into the definition (4.1) of the free energy density, we obtain / = log8 — 
| = 1.32944... as a prediction from the KNS conjecture. 
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